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Abstract 

The value of the internal energy per spin is independent of the strip width for a certain class of 
spin systems on two dimensional infinite strips. It is verified that the Ising model on the kagome 
lattice belongs to this class through an exact transfer-matrix calculation of the internal energy for 
the two smallest widths. More generally, one can suggest an upper bound for the critical coupling 
strength K c {q) for the g-state Potts model from exact calculations of the internal energy for the two 
smallest strip widths. Combining this with the corresponding calculation for the dual lattice and 
using an exact duality relation enables us to conjecture the critical coupling strengths for the three- 
and four-state Potts models on the kagome lattice. The values are K c (q = 3) = 1.056 509 426 929 
and K c (q = 4) = 1.149 360 587 229 2, and the values can, in principle, be obtained to an arbitrary 
precision. We discuss the fact that these values are in the middle of earlier approximate results 
and furthermore differ from earlier conjectures for the exact values. 
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The finite-size scaling technique, known as phenomenological renormalization, has proven 



to be a very reliable method for obtaining critical properties of low- dimensional systems 



Since its beginning, this method has been used to extract thermodynamic quantities in the 



infinite-width limit from transfer- matrix calculations for infinite strips of finite width 
Since a transfer-matrix calculation gives exact results for an infinite strip with width L, 
it can give very precise information on the system in the thermodynamic limit, provided 
that the relevant thermodynamic functions have good convergences with respect to the 
width L. For this reason, corrections to the critical finite-size scaling have been a key issue 

nn 

in this phenomenological renormalization method |4J, |5[. The present investigation makes 
use of the recent progress in computing algorithms, which makes it possible to solve an 
eigenvalue problem for a transfer matrix with a considerable size in an exact or arbitrarily 
precise manner. This means that one may use symbolic algebra systems to solve a given 
transfer matrix in a closed form. Alternatively, if this is not possible, one can do numerical 
calculations with an arbitrary precision, which means that the numerical precision of every 
calculation is free from rounding errors but limited only by the computing memory. This 
makes it possible to take full advantage of the exactness of the transfer-matrix method. In 
the present work we use exact calculations of narrow infinite strips to locate the critical 
point of the g-state Potts models \(^ on the kagome lattice with q = 3 and 4. The kagome 
lattice is one of the simplest two-dimensional (2D) structures belonging to the Archimedean 
lattices and has also drawn practical attention due to distinct structural properties!?, 8]. 
The case of q = 2 on the kagome lattice was solved more than half a century ago [9|], but 
the three- and four-state Potts models have been long-standing open questions in statistical 



physics, and have given rise to, by now, classical 



numerical and approximate determinations 



12 



conjectures 



10 



-[l5|. 



llj , as well number of 



The zero-field g-state Potts model is defined by the following Hamiltonian: 

H=-jY,KSi,Si), 

(ij) 

where each spin S& may take an integer value from to (q — 1), 5 denotes the Kronecker 
delta function, and the sum is over all the nearest-neighbor pairs. We will set the interaction 
strength J as unity throughout this work and identify the inverse temperature /3 with the 



coupling strength K = (3 J. According to the Fortuin-Kasteleyn representation [16|, the 
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partition function corresponding to this Potts Hamiltonian can be written as 

{S} {S} 

with p = 1 — e~ K , where the sum is over all the spin configurations with N c clusters made 
of b connected bonds out of B total bonds inside the system. At the critical point K c the 



partition function for the infinite system (= 
in its .fT-derivatives. The conjecture in Ref. 



Doth length and width infinite) has singularities 



10| states that the critical points can be located 



by solving the following sixth-order polynomial: 

v G + 6v 5 + 9v 4 - 2qv 3 - 12qv 2 - 6q 2 v - q 3 = (2) 



with v = e — 1. As will be described below, our estimates, based on the two thinnest 
infinite strips, are very close to the values predicted by this conjecture. 

To illustrate the transfer-matrix method 17|, we first consider a thin strip of spins with 



size oo x L as shown in Fig. Ufa). Once the transfer matrix is obtained, the free energy per 
spin fi is given in terms of the largest eigenvalue A^ of the matrix as 

- /3f L = L' 1 log \f\ (3) 

and the internal energy per spin is therefore given as 

Furthermore, given the second largest eigenvalue , the inverse correlation length is ob- 

_1 A (0) 

tained as ^ L = log j^fej- m ^ ne context of the arbitrary-precision arithmetic, the differen- 
tiation in Eq. (jlj) may need some care. The derivatives of the eigenvalues can be calculated 
by using the equation A' = Y*T'X, where T' = dT/d/3 is the first-order derivative of the 
matrix T and A is a diagonal matrix with the eigenvalues of T [18l |. The matrices X and 
Y represent the right and left eigenvectors, respectively, which are constructed in such a 
way that Y*X = I is the identity matrix. Here, the asterisk * means the complex con- 
jugate transpose. Suppose that the spins on the strips are described by the g-state Potts 
model with q = 2, which is equivalent to the Ising model with the temperature divided by 
2. It has been shown in Ref. jjj that the correlation length is very well approximated by 
(i oc I near the critical coupling strength K c = f3 c J = log(l + a/2) Ji], 20]. This means 



that close enough to the critical coupling strength one may use the width L as a substi- 
tute for the correlation length and describe the system in terms of this length scale. The 
proportionality coefficient between £^ and L in the limit L — » oo is related to the correlation- 



decay exponent n by conformal invariance 



2 ill. In fact, it has been furthermore found that 



a t-i 



-2L £ 

d/3 V 



P=Pc 



£l {d&fdfi) = const, for every finite L [221 ] . Assuming that £l ~ (0 — 0r) " 



22]. 



near /3 C in the limit of L — > oo, this yields the exact correlation-length exponent v — \ 
Another interesting fact, crucial for the present investigation, is that the internal energy per 
spin [Eq. (j3J)] has at ft — fl c the same value for all the strips irrespective of their widths 



23]. This fact opens up a simple and practical way of locating the critical point of the 2D 



Ising model by calculating the internal energy for the two thinnest strips and then finding 
the coupling strength for which they have the same internal energy. For the square-lattice 
strip in the diagonal direction shown in Fig. Hfb), for example, equating the internal energy 
per spin for L = 2 to that of L = 3, we get 



e 2K (e 2K 



l)(e 



IK 



Qe 2K + 1) 



{e 2K + l)F(K) + (e 2K + l) 2 ^F\K) 



0. 



where F(K) = e — 8e +30e — 8e +1. It is straightforward to see that the nonnegative 
solutions of the equation are e K = 0, 1, oo, and V2±l. Only the latter two are nontrivial and 
give us the exact K c for the 2D ferromagnetic and antiferromagnetic models, respectively. 
The largest eigenvalues for L = 4 and 5 are also available in closed forms and lead to the 
same conclusion. By using the geom etry shown in Fig. Ufa), the exact results e K = 1 + yfq 
have been obtained for q < 5 23j, but there are also cases where the method does not 



apply 23j, |24j . The invariance of the internal energy with respect to strip width has therefore 
been conjectured to be due to certain symmetries in the model 24 . 125 1 . 



We apply the transfer-matrix method to the Potts models on the kagome lattice. First, 
we verify that the known exact solution for the two-state Potts model is reproduced by 
assuming that the internal energy per spin is invariant also in this case. We construct two 
spin blocks for generating the kagome lattice, as illustrated in Figs. [2](a) and^b). Writing 
down the corresponding transfer matrices and denoting their largest eigenvalues as and 
a£°\ respectively, we compute the internal energies per spin as u a = — (3Ai°' ) )~ 1 (9Ai ^/9 / 9) 
for Fig. Eta) and w fe = — (QX^)~ 1 (d\^ /d/3) for Fig. E^b). Indeed, it is readily found that 
u a = ut, = —(7 + 2y / 3)/6 at f3 — ~ log(3 + 2v^3), which is the exact critical point of this 
system jsj. This verifies that the correct critical K c can be obtained from the two thinnest 




FIG. 1: Spin blocks to make spin strips of (a) the square-lattice type, and (b) the double-square- 
lattice type. The periodic boundary condition is imposed in the vertical direction for all the cases, 
so the vertical lengths are regarded as L = 3. Note that the periodic boundary condition may 
introduce double connections in some pairs of spins if L is small. 



A B 




FIG. 2: Spin blocks to make spin strips of the kagome type with (a) L = 1 and (b) L = 2, and 
those of the dice type with (c) L = 1 and (d) L = 2. The dotted lines show the periodic boundary 
condition in the vertical direction. Note that B and C are doubly connected in (a), and such double 
connections are also found in (c). 

strips also for the two-state Potts model on the kagome lattice. Or, in other words, this 
shows that the internal energy per spin is invariant also for the two-state Potts model on 
the kagome lattice. 

In case of the three- and four-state Potts models on the kagome lattice, neither the critical 
values K c are exactly known, nor is it a priori known if the internal energy is invariant. In 
order to generalize the method, we study pairs of strips, one with finite length M and width 
L and the other with 2M and 2L. The two strips in such a pair have the same aspect ratio 
r = M/L where M and L are chosen such that r is a positive number. When r is suitably 
chosen, each such pair will have a single K cvoss (L) for which the internal energy per spin is the 
same. Figures [3](a) and[3](b) show that the crossing point K cross (L) monotonically decreases 
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FIG. 3: Monte Carlo data of the internal energy per spin for q = 3, with aspect ratios (a) r = 2 
and (b) r = 4, respectively. The vertical dotted lines show the critical coupling strength obtained 
in this work, (c) Sectional view at this particular point, where the horizontal dotted line indicates 
the value obtained in this work, u = —1.629 543 706 399 6. The scaling dimension x is expressed 
as (1 — a)/v with the specific-heat exponent a, and therefore x = 4/5 for q = 3. (d) Qualitatively 
the same behavior is observed for q = 4. Here we plot it with r = 4. 

with increasing L. In the limit of L — > oo, the crossing point approaches the true critical 
K c , that is, [K CTOSS (L) — K c ] — >■ + . This implies that K CTOSS (L) will give an upper bound 
for K c for each of the fixed aspect ratios. As long as the aspect ratio does not change any 
essential physics but only the convergence rate toward the bulk criticality [Fig. [3](c)], we can 
suggest that the crossing point K CTOSS (L = 1) for r = oo will either give an upper bound or 
alternatively the exact results: that is, for the case of q = 2 it gives the exact result, whereas 
it gives at least an upper bound for q = 3 and q = 4. An argument can be given in the 
following way: the crossing point would fail to be an upper bound if crossing could be found 
on both sides of the true critical point K c . This actually means that the internal energy per 
spin ul would not be a monotonic function of L at K = K c . We note that the classical Potts 
model on the L x M strip can be mapped to the one- dimensional quantum Potts model of 
size L by putting the strip length M in the imaginary-time direction |26|. One can describe 
the finite-size scaling around the critical point as ul — u^, = L~ x a \(K — K^L 1 ^" , ML~ Z ] 
to the leading order with a two-parameter function a 271 ] . the dynamic critical exponent 
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z — \ |26j, and the scaling dimension x of the energy- density operator [28] . At K = K c , 
the scaling function reduces to a(ML~ l ) = a(r), so we find that ul ~ «oo + a{r)L~ x . The 
geometric factor r is absorbed by the coefficient a which determines the convergence rate. 
It would be plausible to say that a(r) is continuous and nonvanishing for any finite r and 
hence cannot change the sign. The theory of the finite-size scaling therefore tells us that 
ul is a monotonic function of L so that the crossing point for a fixed aspect ratio will exist 
only on one side, which is K > K c in this case. 

In the present investigation we are for practical reasons restricted to q < 4 and L < 2 
since the transfer-matrix size increases as q 3L x q 3L . It is straightforward to write the transfer 
matrices for q = 3 and 4 and solve the eigenvalue problem. By equating u a to Uf, as above, we 
obtain the two values K cvoss (L = 1), which are 1.056 509 426 929 and 1.149 360 587 229 2 
for q = 3 and q — 4, respectively. These values are shown in Table [J together with other 
existing estimates. Note that the arbitrary-precision arithmetic can make our values as 
precise as we want, in principle. As seen in Table [H our values are somewhere in the middle 
of the earlier existing estimates and conjectures, suggesting that they may be the exact 
values. In order to further examine the obtained values, we make use of the fact that there 



exists an exact relation between the critical coupling strengt 
dual [called a dice lattice, compare Figs. [2](c) and[2](d)] [10 



is of the kagome lattice and its 



29J 



(e*--l)(eP--l) = q, (5) 

where K c means the critical coupling strength of the dice lattice. This means that the upper 
bound K cvoss (L = 1) obtained for the dual lattice can be turned into a lower bound for K c 
of the kagome lattice. Repeating the calculation for K cmss (L = 1) with the two thinnest 
strips, given in Figs. [2(c) and|2J(d), gives K CTOSS (L = 1) = 0.955 080 568 397 4 as an upper 
bound for K c , which through Eq. (J5j) gives the lower bound K c = 1.056 509 426 929 0. 
This is, to all the 14 decimal places, identical to the upper bound obtained directly for the 
kagome lattice. The most reasonable conclusion is that the calculation gives the exact value 
and that, just as for the two-state Potts model on the kagome lattice, the internal energy is 
independent of the strip width at the critical temperature and that, furthermore, the same 
is true for the q = 4 case [Fig. [3(d)]. A complete analytic argument is called for, and a 
simple way to test this conjecture would be to solve the transfer matrix with L = 3. 

As seen from Table [J the situation for the three-state Potts model is as follows: both 
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TABLE I: Critical thresholds of the (/-state Potts models on the kagome lattice in terms of p 

l-e~ K . 

Reference <7 = 1 q = 2 q = 3 « = 4 



Exact [9] 1 - 1/V3 + 2^/3 



Conjecture [1Q] 0.524 429 71 1 - 1/v 3 + 2^3 0.652 327 40 0.683 127 34 



Conjecture [11] 0.522 372 07 1 - 1/V3 + 2^3 0.653 932 82 0.685 967 83 

Series [12] 0.652 350(5) 0.683 15(5) 
Monte Carlo [30] 0.524 405 3(3) 

Monte Carlo [13] 0.606 62(8) 0.652 32(7) 0.683 17(2) 

Subnet [15] 0.524 404 978(5) 0.606 680 106 83(15) 0.652 350 2(4) 0.683 163(5) 



This work 1 - 1/V3 + 2^/3 0.652 332 747 264 01 0.683 160 704 284 84 



the two earlier conjectured exact values can be ruled out, although the conjecture by Wu in 
Ref. [10] is very close to the value in this work. Our conjectured exact value is somewhat 
surprisingly outside the bounds of the value estimated from series expansion in Ref. [12[ and 
the subnet estimate in Ref. 15| . It agrees well with and is inside the bounds of the Monte 



Carlo estimate in Ref. 



131 ] . For the four-state Potts model, the situation is somewhat differ- 



ent: again, the earlier conjectured exact values can be ruled out. However, our conjectured 
exact value is inside the bounds of all the other estimates. 

The conjecture by Wu in Ref. [h]] gives the critical coupling strengths as solutions of the 
sixth-order polynomial given in Eq. ([2]). It is important to note that we have also given 
our values as solutions of certain polynomial equations since we are dealing with transfer 
matrices. Although we have not factorized the full polynomials yet, one may ask if such 
sixth-order pol yno mials as Wu has derived can be eventually factored out. To answer this, 



we follow Ref. 



30j and try to determine polynomials in the variable v = e K — 1 which have 



roots at the exact critical values. Even if we work with numeric values, instead of symbolic 
manipulations, this method makes it possible to find such a polynomial. For example, in 



case of q — 2, one recovers the compact analytic expression v = y 3 + 2-^/3 — 1 by solving 
the obtained polynomial equation. Based on the conjecture by Wu, we try to find the value 
for the q = 3 case as the solution of the sixth-order polynomial Ylt=o c n v n = with integer- 
valued coefficients. We furthermore assume that Cq = 1, and let C5 and C4 vary from —25 



S 



to 25, while the other four coefficients may take values from — 10 2 to 10 2 . Substituting our 
value v = 1.876 313 463 895 for q = 3, the best polynomial is found to be 

v 6 - 6t> 5 + 22t> 4 - 79v 3 + 99v 2 + 28v - 56 = 0, 

yielding a solution v poly = 1.876 313 463 898. Even if the discrepancy between our con- 
jectured exact value and the solution of the polynomial is tiny, it is still significant, which 
means that there is no such polynomial within the range of coefficients tested. This might 
suggest that the solution cannot be obtained from a simple sixth-order polynomial as was 
assumed in the conjecture by Wu. 

We conclude this work with a brief sideline: the Fortuin-Kasteleyn representation [Eq. (pQ)] 
for q = 1 recovers the bond-percolation problem on the kagome lattice. The method in the 
present paper cannot be directly used in this case, since the internal energy becomes a 
constant independent of the coupling strength. Without the knowledge of the exact in- 
dependence of K c , one can only interpolate it from the other estimates. For the case of 
the square lattice, which has coordination number 4, as does the kagome lattice, we have a 
general expression of the critical point as t> sq = yfq [6|. We assume that v(q) of the kagome 
lattice can be expanded in series of this variable: v(q) = a(y/q) 3 + b(y/q) 2 + Cy/q, where we 
further note that v(0) = is an exact limit. Finding the three parameters a, b, and c by 
substituting the conjectured values for v(2), v(3) and v (4), we can interpolate the value at 
q — 1 and obtain p « 0.52433. Compared to the numerical estimate shown in Table HI the 
fractional error amounts to be about 140 parts per million. 

In summary, we have conjectured the exact values of critical temperatures for the three- 
state and four-state Potts models on the kagome lattice by using exact transfer-matrix 
calculations on thin infinite strips. This suggests that the internal energy can provide a 
sharper condition for criticality than lattice symmetries considered in the earlier conjecture 
by Wu. It has also been noted that, for the three-state Potts model on the kagome lattice, 



the series expansion in Ref. [12| does not contain our result within its bounds. The method 
devised to obtain the results is based on exact solutions of the two thinnest infinite strips. 
These solutions have been obtained by taking full advantage of computational symbolic 
algebra systems. Since the method itself appears to be quite general, it may possibly be 
used to solve other problems. 



9 



Acknowledgments 



S.K.B. and P.M. acknowledge the support from the Swedish Research Council with Grant 
No. 621-2008-4449. H.M. thanks E. Lundh for access to computational resources. B. J.K. was 
supported by Basic Science Research Program through the National Research Foundation of 
Korea (NRF) funded by the Ministry of Education, Science and Technology (2010-0008758). 
This research was conducted using the resources of High Performance Computing Center 
North (HPC2N). 



[1] P. Nightingale, J. Appl. Phys 53, 7927 (1982). 

[2] M. E. Barber, in C. Domb and J. L. Lebowitz, eds., Phase Transitions and Critical Phenomena 

(Academic Press, London, 1983), p. 145. 

[3] T. W. Burkhardt and J. M. J. van Leeuwen, Real-Space Renormalization (Springer, Berlin, 

1982). 

[4] J. M. Luck, Phys. Rev. B 31, 3069 (1985). 

[5] S. Singh and R. K. Pathria, Phys. Rev. Lett. 55, 347 (1985). 

[6] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 

[7] A. Souslov, A. J. Liu, and T. C. Lubensky, Phys. Rev. Lett. 103, 205503 (2009). 

[8] Q. Chen, S. C. Bae, and S. Granick, Nature (London) 469, 381 (2011). 

[9] K. Kano and S. Naya, Prog. Theor. Phys. 10, 158 (1953). 

[10] F. Y. Wu, J. Phys. C 12, L645 (1979). 

[11] C. Tsallis, J. Phys. C 15, L757 (1982). 

[12] I. Jensen, A. J. Guttmann, and I. G. Enting, J. Phys. A 30, 8067 (1997). 

[13] J.-A. Chen, C.-K. Hu, and F. Y. Wu, J. Phys. A 31, 7855 (1998). 

[14] F. Y. Wu, Phys. Rev. E 81, 061110 (2010). 

[15] C. Ding, Z. Fu, W. Guo, and F. Y. Wu, Phys. Rev. E 81, 061111 (2010). 

[16] C. M. Fortuin and P. W. Kasteleyn, Physica 57, 536 (1972). 

[17] J. J. Binney, N. J. Dowrick, A. J. Fisher, and M. E. J. Newman, The Theory of Critical 

Phenomena (Oxford University Press, New York, 1992). 

[18] N. van der Aa, Perturbation theory of eigenvalue problems (2005), seminar in Centre 



10 



for Analysis, Scientific computing and Applications, Technische Universiteit Eindhoven, 
http : //www . win . tue . nl/ casa/meetings/ seminar/previous/abstract 05 10 19 . html, 

[19] H. A. Kramers and G. H. Wannier, Phys. Rev. 60, 252 (1941). 

[20] L. Onsager, Phys. Rev. 65, 117 (1944). 

[21] J. L. Cardy, J. Phys. A 17, L385 (1984). 

[22] M. A. Yurishchev, Nucl. Phys. B Proc. Suppl. 83-84, 727 (2000). 

[23] J. Wosiek, Phys. Rev. B 49, 15023 (1994). 

[24] I. Souza, C. Lobo, and J. Seixas, Phys. Rev. B 55, 6356 (1997). 

[25] M. A. Yurishchev, Nucl. Phys. B Proc. Suppl. 106, 917 (2002). 

[26] T. Senthil and S. N. Majumdar, Phys. Rev. Lett. 76, 3001 (1996). 

[27] H. Rieger and A. P. Young, Phys. Rev. Lett. 72, 4141 (1994). 

[28] M. Weigel and W. Janke, Phys. Rev. B 62, 6343 (2000). 

[29] F. Y. Wu, J. Math. Phys 18, 611 (1977). 

[30] R. M. Ziff and P. N. Suding, J. Phys. A 30, 5351 (1997). 



11 



